Controllable magnetic doping of the surface state of a topological insulator 
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A combined experimental and theoretical study of doping individual Fe atoms into Bi2Se3 is pre- 
sented. It is shown through a scanning tunneling microscopy study that single Fe atoms initially 
located at hollow sites on top of the surface (adatoms) can be incorporated into subsurface layers 
by thermally-activated diffusion. Angle-resolved photoemission spectroscopy in combination with 
ab-initio calculations suggest that the doping behavior changes from electron donation for the Fe 
adatom to neutral or electron acceptance for Fe incorporated into substitutional Bi sites. According 
to first principles calculations within density functional theory, these Fe substitutional impurities 
retain a large magnetic moment thus presenting an alternative scheme for magnetically doping the 
topological surface state. For both types of Fe doping, we see no indication of a gap at the Dirac 
point. 



PACS numbers: 68.37.Ef, 71.15.Mb, 73.20.At, 79.60.-i 

Three-dimensional topological insulators (3D TIs) 
have a gapped bulk band structure but a spin-orbit in- 
teraction induced band inversion leads to the existence 
of a metallic topological surface state (TSS). The TSS 
can exhibit a Dirac-like dispersion which is protected by 
time-reversal symmetry, resulting in unique properties 
such as spin-momentum locking and prohibited backscat- 
tering [TJ [2] - Nevertheless, defect scattering is possible 
for a two-dimensional TSS and scattering by magnetic 
impurities is of particular interest, as these can break 
time-reversal symmetry in the system. This may affect 
the topological protection and could potentially open a 
gap in the topological protected band connection at the 
surface [3H6]. 

Two experimental approaches have been utilized to 
introduce magnetic impurities in 3D TIs, namely bulk 
magnetic doping in crystals or thin films [5| f7HTQ] . and 
in-situ surface doping via adsorption of impurities on 
pristine topological insulators [6j [TTJ [12] . Surface dop- 
ing can be easier to control experimentally and has the 
additional advantage that the bulk material remains a 
TI, something that is not necessarily the case for too 
high concentrations of bulk dopants [7] . Angle-resolved 
photoemission spectroscopy (ARPES) measurements of 
magnetically-doped Bi-chalcogenide samples prepared in 
both manners exhibit a feature near the DP which resem- 
bles a gap [5, 6, 8 . For the surface-doped case, however, 
it turned out that similar features can be generated by 
doping with non-magnetic impurities [13 . Indeed, an 



apparent gap opening can be caused via the adsorbate- 
induced strong band bending near the surface p~4] [15] . 
More recently, gaps in bulk magnetically doped thin 
Bi2Se3 films have been observed [8 resulting in a spin 
reorientation of the TSS. Finally, it has been shown that 
potential disorder induced by charged bulk defects hin- 
ders a clear analysis of faint effects due to the breaking 
of time-reversal symmetry on the scattering of the TSS 
electrons [16 . Therefore, a method where neutral mag- 
netic impurities can be introduced into the surface and 
magnetically interact with the TSS, while only weakly 
perturbing the band structure, is highly desirable. 

In this letter, we demonstrate a way to control the 
charge state of a magnetic impurity such that this is 
achieved. Scanning tunneling microscopy (STM) exper- 
iments illustrate that Fe atoms that are adsorbed on 
the surface (adatoms) diffuse into subsurface sites of 
Bi2Se3 upon thermal annealing. Such bulk doping via 
the surface can be realized for a large range of coverages. 
Complementary ARPES studies illustrate that annealing 
reduces the downwards band-bending, initially produced 
by Fe adatoms, and restores the DP near its initial loca- 
tion before Fe adsorption. These effects are seen regard- 
less of the initial location of the DP, namely for both 
pristine (n-doped) and Ca-doped (nearly intrinsic) sub- 
strates. The difference in doping behavior of Fe adatoms 
compared to incorporated subsurface Fe atoms at dif- 
ferent lattice sites is further studied by first principles 
calculations within density functional theory. The cal- 
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FIG. 1: Constant current STM images of differently prepared samples of Fe on Bi2Se3 (It = 50pA, Vt> = 0.25V (a-d), Vb = IV 

(e) .) (a) Cold deposited Fe (~ 1% of a ML) on pristine Bi2Se3. Inset: magnified view showing 8 individual Fe adatoms (lOnm 
x lOnm, Vb = —0.3V, It = InA). (b), (c) Samples prepared as in (a) annealed to Ta ~ 260 K and 370 K, respectively, (d) Cold 
deposited Fe (0.5% ML) on Ca-doped Bi2Se3. (e) Sample prepared as in (d) but with 6% ML Fe and annealed to Ta ~ 370 K. 

(f) Number of Fe adatoms and subsurface Fe atoms relative to initial Fe adatom number as a function of Ta. Emergent clusters 
were not counted in the statistics. 



dilations reveal that the doping character changes from 
electron donor for Fe adatoms towards neutral for Fe in- 
corporated into Bi substitutional sites consistent with the 
ARPES results. Such incorporated Fe atoms theoreti- 
cally retain a large magnetic moment which is a prerequi- 
site for modifying the TSS by breaking the time-reversal 
symmetry locally. 

STM experiments were carried out in a multichamber 
UHV-system with a base pressure below 1 x 10 -10 mbar 
utilizing a home-built variable temperature scanning tun- 
neling microscope (VT-STM) [17] • Annealed tungsten 
tips were used for all experiments and both tip and sam- 
ple were cooled to T = 30 K. STM topography images 
were taken in constant current mode at a tunneling cur- 
rent I t with the bias voltage Vb applied to the sample. 
ARPES measurements were performed at the SGM-3 
beamline of the ASTRID synchrotron radiation facility 
[T8] at a sample temperature of T = 60 K. The angular 
and combined energy resolutions were 0.13° and better 
than 15 meV, respectively. Spectra were obtained along 
the K — V — K azimuthal direction with a photon energy 
of hv = 19.2 eV. Pristine and Ca-doped Bi2Se3 single 
crystals were grown as described in refs. [l4j[T9]. In both 
experimental setups the Bi2Se3 samples were cleaved in 
situ at room temperature (RT) and immediately cooled 
to temperatures of T < 150 K. Fe was originally de- 



posited onto the cold surface by e-beam evaporation re- 
sulting in a distribution of single Fe adatoms [20] . 

To achieve a theoretical description of the Fe/Bi2Se3 
system, we have performed density functional theory cal- 
culations with the VASP code [2TJ [22], where we em- 
ployed the GGA-PBE approximation to the exchange 
and correlation functional 23 and used plane wave basis 
sets with a kinetic energy cutoff of 20 Ry in combina- 
tion with projector augmented- wave (PAW) pseudopo- 
tentials [24j [25]. We used a (3 x 3) supercell with two 
quintuple slabs which have been relaxed until the max- 
imum force component of every atom was lower than 
0.02 eV/A. A (6 x 6) Appoint grid has been employed 
to obtain the precise density of states (DOS). 

Fig. [I] shows a sequence of STM images where Fe was 
initially deposited onto the cold surface (T < 150 K) of 
pristine (a-c) and Ca-doped (d,e) Bi2Se3. After depo- 
sition, the sample was annealed to a given temperature 
Ta as indicated in the image. For each subsequent an- 
nealing temperature, the sample was re-cleaved and the 
process repeated. After cold deposition (Fig. [IJa) and 
inset), STM reveals that the surface is covered primar- 
ily by individual Fe adatoms which are strongly relaxed 
into either of the two hollow sites of the hexagonal lat- 
tice of Se surface atoms like seen in ref. [20]. After an- 
nealing to T A « 260 K (Fig. [TJb)), a new type of de- 
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FIG. 2: (a), (b) Atomically resolved STM images of two 
different subsurface Fe defects [it = 50pA, Vb = IV (a), V\> = 
0.25V (b)]. (c)-(h) Bias dependent constant current images 
of a subsurface Fe atom (red arrows), 3 Fe adatoms and a Se 
vacancy at Vb as indicated, (i t = lOOpA, Az = 35pm). 



feet appears, which is imaged as a triangular-shaped de- 
pression at Vb = 0.25 V. Simultaneously, a decrease in 
the overall density of Fe adatoms and the formation of 
clusters are observed (larger bright spots in Fig. [TJb)) . 
After annealing to T A « 370 K (Fig. [jjc)) the den- 
sity of the triangular-shaped depressions increases signif- 
icantly accompanied by an increase in the size of the clus- 
ters, at the expense of the number of Fe adatoms. The 
same qualitative behavior is found for Ca-doped sam- 
ples (Figs. [TJd,e)). However, we focus our quantitative 
analysis on Fe on the pristine Bi2Se3 surface since the 
triangular-shaped depressions look similar to those of Ca 
dopants [26] . 

In Fig.[TJf) the number of Fe adatoms and triangular- 
shaped depressions relative to the initial number of de- 
posited Fe adatoms, as a function of Ta, is depicted. 
The abrupt decrease in the number of Fe adatoms at 
Ta ~ 275 K coincides with the emergence of the de- 
pression defects. This strongly indicates that the ob- 
served triangular-shaped depressions are the signature 
of Fe atoms incorporated into a subsurface layer of the 
Bi2Se3 substrate, upon thermal activation. In order to 
further exclude that the triangular-shaped depressions 
originate from thermally-induced defects of the host ma- 
terial, e.g. from Bi vacancies, we did the following con- 
trol experiment: Fe was initially deposited at RT and 
subsequently cooled down for STM measurements (not 
shown) revealing a similar number of triangular-shaped 
depressions as in the case where cold deposition and sub- 
sequent RT annealing were performed. The density of 
triangular-shaped depressions correlates with the amount 
of deposited Fe as shown for a higher initial coverage of 
6% ML in Fig. [TJe). Moreover, the clean Bi2Se3 surface 
without the presence of Fe does not exhibit such defects 
regardless of applied temperature. 

Closer analysis of the STM topographs provides a 
deeper insight into the subsurface doping process. For 



Ta ~ 275 K (Fig. [TJb)), only one type of triangular- 
shaped depressions is found. Samples that were annealed 
at higher temperatures exhibit triangular-shaped depres- 
sions with at least two different appearances in STM 
topographs (Fig. [ljc)) . Figs. j2^a,b) exemplarily show 
high resolution STM images of two such defects and 
Figs. |2|c-h) show bias-dependent STM images of one of 
those defects. Defects with a strikingly similar STM ap- 
pearance have been recently observed for Fe bulk-doped 
Bi2Te3 [27], Fe adatoms on Bi2Te3 after annealing [28] 
and for Fe bulk-doped Bi2Se3 [29]. Such defects were 
attributed to substitutional Fe atoms within Bi sites in 
the first quintuple layer (QL), i.e. in the 2 nd and 4 th 
subsurface layers. We thus conclude the observed tri- 
angular depressions are substitutional Fe atoms residing 
in Bi sites. Moreover, as pointed out in Ref. [29], Fe 
atoms substituted into Bi sites in deeper QL are too far 
below the surface to be imaged by STM. This might ex- 
plain part of the discrepancy between the deposited Fe 
and the number of Fe substitutional defects counted af- 
ter annealing (Fig.[TJf)). Besides, emerging clusters upon 
annealing, which may contain both Bi and Fe, may also 
account for this discrepancy. 

To further investigate the impact of Fe doping on the 
TSS, ARPES measurements were performed. In Fig. [3] 
two series of ARPES measurements are shown for pris- 
tine (a-e) and Ca-doped (f-j) Bi2Se3 substrates, at dif- 
ferent stages of the preparation procedure. After in- situ 
cleaving (a,f), subsequent iterations of cold Fe deposi- 
tion (T « 150 K) and annealing to Ta ~ 370 K were 
performed with no cleaving in between. The amount 
of subsequently deposited Fe was increased for each it- 
eration (1%, 2.5%, 5%, 10%). The ARPES data af- 
ter each cold deposition but before subsequent annealing 
(Fig. [3] (b,d,g,i)) reveals a downward shift of the DP as 
compared to its energetic position before Fe deposition. 
We can thus conclude that Fe adatoms are donors in- 
ducing downwards band bending. Such band bending 
leads to the well known formation of quantum well states 
(QWS) [6] [LQ USUI EGS [30] in both the conduction and 
the valence bands [13 , where the conduction band QWS 
shows a clear Rashba splitting. 

Surprisingly, after annealing to temperatures coincid- 
ing with STM observations of Fe diffusion into the sub- 
surface region (Fig. T|2), the DP shifts back toward its 
initial position and the surface's electronic band structure 
is reverted close to its initial state before Fe deposition 
(Fig. [3] (c,e,h,j)). Strikingly, after annealing, the QWS 
are nearly destroyed. In order to further analyse this be- 
havior, the binding energy of the DP after each prepara- 
tion step of the subsequent iterations has been extracted 
from the ARPES spectra and is plotted in Fig. [3jk) for 
the pristine and two Ca-doped substrates with different 
Ca concentrations. It reveals the afore mentioned oscil- 
lation of the DP binding energy after each preparation 
step, i.e. an increase in the DP binging energy for cold 
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FIG. 3: Two series of ARPES measurements for pristine (a-e) and Ca-doped (f-j) Bi2Se3 samples after different Fe deposition 
and annealing steps, (a), (f) Initial spectra after in-situ cleaving. Subsequent cycles of cold deposition of Fe (T « 150K, 
% ML Fe as indicated and annealing at Ta ~ 370K followed. The DP extracted from the spectra is marked by a green and 
red horizontal line, respectively, (k) The resulting binding energy of the DP is plotted in dependence of the preparation step 
number for a pristine (circles) and two differently Ca-doped samples (up and down triangles). Preparation step number 
indicates the freshly cleaved substrate without Fe, preparation step 1 (1%), 3 (2.5%), 5 (5%), 7 (10%), and 2, 4, 6, 8 subsequent 
annealing steps. 



deposition of Fe, and a decrease after annealing the Fe 
covered sample. The gradual shifting of the DP toward 
higher binding energies is most likely the result of the 
well known aging effect of the Bi2Se3 surface [30 . In 
order to exclude the influence of rest gas physisorption 
at low temperature that would also lead to a downward 
band bending, the following test has been performed. 
The very same procedure of Fe deposition has been fol- 
lowed but keeping the shutter of the Fe evaporator closed 
and no appreciable change in the band bending has been 
observed after this simulated deposition. Our result is 
contrary to a recent ARPES study that concludes sur- 
face electron doping for Fe deposition at RT and hole 
doping for cold deposited Fe [12 . Taking into account 
the STM-observed thermally-induced subsurface diffu- 
sion of Fe adatoms (Fig. [I] (f)), our sequence of ARPES 
experiments indicates, that substitutional Fe atoms are 
either neutral dopants or electron acceptors, while the Fe 
adatoms that are prepared by cold deposition are elec- 
tron donors. In all of these cases the DP stays intact 
without the opening of a gap. 

In order to determine the electronic and magnetic 
properties of the incorporated Fe atoms, we performed 
ab initio calculations. As shown in the ball and stick 
models in Fig.[4](a-f), we considered pure Bi2Se3, as well 
as Fe on the two different hollow sites fee (a) and hep 
(b), Fe substituting Bi atoms in the 2 nd (c) and 4 th (d) 
subsurface layer, and Fe in two different interstitial sites 
(e,f) of the Bi2Se3 lattice. The calculated total density 
of states (DOS) for these different cases is illustrated in 
Fig. [4] (g). For pure Bi2Se3 (grey shaded), the onset of 



the valence band emerges at energies below the Fermi en- 
ergy Ey and the onset of the conduction band DOS for 
E — Ey > 0.75 eV. The minimum and gradual increase in 
the DOS within the bulk band gap (E - E F < 0.75 eV) 
are due to the DP of the TSS. Compared to the DOS of 
the pure Bi2Se3, the Fe adatoms (green, yellow) or the Fe 
subsurface interstitials (cyan and purple) cause a shift of 
the bulk band gap and DP to a higher binding energy by 
AE w 0.5 eV. This reproduces the electron doping be- 
havior of the Fe adatoms, as observed by ARPES. More- 
over, Fe atoms in Bi substitutional sites (red, blue) do 
not shift the bulk band gap and DP with respect to the 
pure Bi2Se3, indicating neutral doping for these impuri- 
ties. Note here that the Fermi energy is always below 
the onset of the conduction band within the calculations, 
while in the experiment it is slightly shifted into the con- 
duction band due to Se vacancy doping and aging effects. 
For a strong shift of the Fermi energy, the Fe atoms on 
Bi substitutional sites might become electron acceptors. 
These conclusions are consistent with the above interpre- 
tation of the STM and ARPES data, confirming that the 
thermally-induced subsurface Fe doping results in a pref- 
erential occupation of Bi substitutional sites, as already 
suggested by STM. 

Finally, we calculated the magnetic spin moments of 
the Fe atoms in the different lattice sites as indicated in 
the legend of Fig. [4](g). For all defects, we find consider- 
able magnetic moments 2.7/i# < ^Fe < 3.5/i# m the Fe 
d shell. mp e is largest for Fe substitutional defects, where 
Fe is nominally in a 3 + oxidation state. In contrast, Fe 
adatoms are clearly in a lower oxidation state as no Bi 3+ 
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FIG. 4: (a-f ) Side views of ball and stick models of the differ- 
ent Fe doped Bi2Se3 systems considered in the ab-initio calcu- 
lations: (a,b) Fe on the two different hollow sites fee (Fe Tl) 
and hep (Fe_T2), (c,d) Fe substituting Bi atoms in the 2 nd 
(Fe_2Bi) and 4 th (Fe_4Bi) subsurface layer, and (e,f) Fe in two 
different interstitial sites (Fe ll, Fe_I2) of the Bi2Se3 lattice. 
Fe is shown in black, Bi grey, and Se orange, (g) Calculated 
total DOS for the different considered sites. The resulting 
energy shift AE of the DOS due to the doping is indicated. 
The Fe 3d magnetic moments and the total magnetic mo- 
ments (mFe, Tritot) for Fe at these different sites are given in 
the legend in units of \±b- 

is missing. Thus, the reduced magnetic moments of Fe 
adatoms as compared to substitutional Fe atoms as well 
as the stronger tendency of Fe adatoms to donate further 
charge to the Bi2Se3 conduction band is well understand- 
able. Fe substitutional atoms thus not only minimize the 
charge doping but also maximize the magnetic coupling 
to the Bi2Se3 host. 

For a magnetism-induced modification of the TSS by 
Fe, the coupling of the Fe 3d magnetic moments to the 
Bi2Se3 host is decisive. It is therefore illustrative to an- 
alyze also the total spin magnetic moments m to t of the 
supercell. While the total magnetic moments m to t and 
mp e are rather similar in the adatom case, m to t exceeds 
m Fe by ~ 1-5/ie in the case of Fe substitutionals. In 
the latter case a sizable magnetic moment is induced in 
the vicinity of Fe which indicates much stronger mag- 
netic coupling between substitutional Fe atoms and the 
Bi2Se3 host than for adatoms. 



In summary, it was demonstrated that, via low temper- 
ature deposition and moderate annealing, neutral mag- 
netic moments can be introduced into the topmost lay- 
ers of a prototypical TI, thus preserving the bulk elec- 
tronic properties of the TI without causing significant 
band bending effects at the surface. In the investigated 
case of Fe on Bi2Se3, Fe acts as an electron donor when 
adsorbed on top of the surface at low temperature, but 
becomes neutral or an acceptor when substituted into a 
subsurface Bi site by thermally- activated diffusion. This 
allows for a study of the time-reversal symmetry break- 
ing driven scattering of the TSS electrons by magnetic 
impurities. 
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